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ABSTRACT 

The spectral energy distributions (SEDs) of some T Tauri stars display a deficit of 
near-IR flux that could be a consequence of an embedded Jupiter-mass planet partially 
clearing an inner hole in the circumstellar disc. Here, we use two-dimensional numerical 
simulations of the planet-disc interaction, in concert with simple models for the dust 
dynamics, to quantify how a planet influences the dust at different radii within the 
disc. We show that pressure gradients at the outer edge of the gap cleared by the planet 
act as a filter - letting particles smaller than a critical size through to the inner disc 
while holding back larger particles in the outer disc. The critical particle size depends 
upon the disc properties, but is typically of the order of 10 microns. This filtration 
process will lead to discontinuous grain populations across the planet's orbital radius, 
with small grains in the inner disc and an outer population of larger grains. We show 
that this type of dust population is qualitatively consistent with SED modelling of 
systems that have optically thin inner holes in their circumstellar discs. This process 
can also produce a very large gas-to-dust ratio in the inner disc, potentially explaining 
those systems with optically thin inner cavities that still have relatively high accretion 
rates. 

Key words: solar system: formation — planets and satellites: formation — planetary 
systems: formation 



1 INTRODUCTION 



A fraction of T Tauri stars, including TW Hva (ICalvet et al. 


2002 


: lUchida et alJ l2004t). GM Aur dMarsh & Mahonev 


1992 


:lKoerner. Sareent & Beckwithlll993l: iRice et alJl2003l). 


CoKu Tau/4 jForrest et ail 12004 lOuillen et alJ 


2004; 


iD'Alessio et alJ l2005|) and DM Tau (ICalvet et al.l 


2005). 



present spectral energy distributions (SEDs) that combine 
little excess flux (above the photospheric level) at near- 
IR wavelengths with robust excesses at A > lO/^m. A 
similar effe ct is seen in at least one Herbig Ae/B e star, 
HD100546 teouwman et al.ll2003t Iciradv et al.ll2005D . Since 



circumstellar discs are generally most optically thick in 
the inner, hottest regions, this observation is noteworthy 
and indicative of some process that can create an opti- 
cally thin inner cavity within the disc. One exciting pos- 
sibility is that the cavity arises from the tidal truncation of 
the inner disc by a massive planet dGoldreich fc Tremaind 
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ll98Ct iLin fc Papaloizoulll986ft . in which case the SED ob- 
servations provide unique constraints on the time scale 
and frequency of giant planet formation. However, non- 
planctary explanations are also possible. Dust can grow 
rapidly to form solid bodies that contribute negligible opac- 
ity in any observed waveband, while larger bodies may 
themselves suffer destructive collisions or rapid orbital de- 
cay due to aerodyna mic drag against the gaseous disc 
JWeidenschillindll977D . Photoevaporation of the outer disc 
by ionizing stellar radiation can also starve the inner disc 
of gas, yielding a short inner hole phase even in t he ab- 
sence of a planet (lAlexander. Clarke fc Pringldl200a) . Pro- 
toplane tary disc models that include s ome or all of these 
effect s l|Taketichj - J^AxtYjnowicj |200l|; |Kejrvmi_fe_13KHnlev, 
| 2004t iDullemond fc Dominikll2005l iTakeuchi. Clarke fc Lin 
2005) admit the formation of dusty rings and inner holes, 
even in the absence of planets and gas. As improved SEDs 
and spectra from Spitzer become available, more detailed 
theoretical study of both classes of model (pure disc and 
disc plus embedded planet) is warranted in order to deter- 
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mine whether these models can be distinguished without 
direct imaging data. 

A sufficiently massive planet (or brown dwarf) can 
probably create an inner cavity in which both the gas 
and dust surface densities are extremely low. However, of 
the four aforementioned T Tauri stars - around which 
inner disc 'holes' are inferred from the lack of near-IR 
emission - only CoKu Tau/4 appears to be entirely de- 
void of circumstellar material within 1 AU of the star. 
In the case o f TW Hya there is b oth ongoing low-level 
gas accretion jMuzerolle et al J 1200(f) . and evidence (from 
the existence of a silicate emission feature at lOpim) for 
a sub-lunar mass of micron sized dust at small radii 
JCalvet et alJl2002l) . Both GM Aur and DM Tau are Clas- 
sical T Tauri stars, with estimated accretion rates on 
to the stellar photosphere variousl y estimated at M ~ 
10" 9 Mcpyr" 1 to 10~ 8 Mr^vr" 1 (iGullbring et al] Il99 " 



Hartmann et alJll998l Iwhite fc Ghezll200ll) . Unfortunately 



Naiita. Carr fc Mathieul ]2003l) failed to detect CO 



emission 

from an inner disc around GM Aur, so there are no obser- 
vational clues as to the gas surface density that must ac- 
company this accretion. However, for these accretion rates 
protoplanetary disc models pre dict gas surface densities at 
r ?a 0.1 AU of 200-2000 gem" 2 jBell et alJll997h . Although 
disc models are subject to large uncertainties arising from 
the unknown efficiency of angular momentum transport 
within protoplanetary discs, the robust accretion signatures 
seen for these stars imply that an inner gas disc must be 
present. 

Within a planetary interpretation of the ob- 
served SEDs, the presence of inner gas discs and 
ongoing stellar accretion is not in itself su rp rising . 
Numerical simulations iArt Ymqwic z fc Lubowl Il996t 
iLubow. Seibert fc Artvmowiczl ll999fT "show that the tidal 
barrier erected by a massive planet is quite leaky. For 
Jupiter mass planets gas can readily overflow the gap, 
and while much of this gas is accreted by the planet in 
excess of 10% of the accretion rate at large radius can 
reple nish the inner disc and ul timately accrete on to the 
star jLubow fc D'Angelol EoOrj . This is consistent with 
GM Aur being a Classical T Tauri star, but the lack of dust 
opacity then requires that the planet needs to be relatively 
more effective at holding back dust - i.e. it needs to filter 
the gas as it passes the tidal barrier. In this paper, we show 
that the required filtration may occur as a consequence of 
the sharp pressure gradient that exists near the outer edge 
of the gap. The pressure gradient leads to super-Keplerian 
rotation of the gas, which in turn allows standard aero- 
dynamic coupling to solid bo dies JWeidenschillind Il977t 
iPaardekooper fc Mellenial l2006f) to inject angular momen- 
tum and hold back the inward flow of all but the smallest 
solid particles. We estimate the efficiency of this process, 
and argue th at it allows relative ly low mass planets - which 
are observed iMarcv et al.ll2005h to be much more common 
than those planets with M p > Mj that would start to cut 
off accretion entirely - to nevertheless create dust-poor 
inner cavities. Within this scenario, the variation in the 
observed properties of T Tauri stars with weak near-IR 
emission could reflect different mass embedded planets - 
relatively low in the case of GM Aur and DM Tau, and 
substantially higher in the case of CoKu Tau/4. 



2 GAS-DUST DYNAMICS 

The disc gas and the embedded dust particles in teract via 
a drag force JWhipplelll972l:IWeidenschilrindll977h that is a 
consequence of the different gas and dust orbital velocities. 
The Keplerian velocity, Vk , at a radius, r, around a star of 
mass M* is 



v K 



GM, 



(1) 



If P is the gas pressure and if p is the gas density, the az- 
imuthal component of the gas velocity, v^, in centrifugal 
equilibrium differs from the Keplerian velocity due to the 
influence of the gas pressure gradient and is given by 



JL 
r 



+ 



ldP 
p dr 



(2) 



The dust particles, on the other hand, are not influenced by 
pressure forces and, in centrifugal equilibrium, orbit at the 
Keplerian velocity, Vk- 

The effect of the drag force depends on the size of the 
dust particles relative to the mean free path of the gas 
molecules. If we assume that the gas is primarily molcec- 
ular hydrogen, the mean free path, A, is 



A 



m,H 2 
pA 



4 x 10" 



(3) 



where mn 2 is the mass of the hydrogen molecule and A is 



nai x 7 x lO'^cm 2 ). The drag force, 



(4) 



its cross-section (A ■ 
Fu is then 

1 ^ 2 2- 

Frj = — — Cuita pu u, 

where u is the relative velocity between the gas and dust 
particles, u = |u|, u = u/u, a is the mean radius of the dust 
grains and Cd is the drag coefficient, given by 

3 u 



a < 9A/4 



C D = i 



24iC 1 



247?; 



0.44 



Re < 1 

1 < R e < 800 
Re > 800. 



(•») 



In equation 0, R e is the Reynolds n umber whi c h can be 
shown to be given by (see for example iRice et al] (120041) ') 



Re 



(6) 



In general, the gas pressure in protoplanetary discs 
decreases with increasing radius. The pressure gradient is 
therefore negative and the gas velocity is sub-Keplerian. The 
drag force then acts to slow the dust particles down caus- 
ing them to spiral in towards the central star. One can then 
determine the radial drift velocity of dust particles by self- 
consistently solving the radial and azimut hal components 
of the momentum equation as outlined in IWeidenschillind 
il977l) . Figure Q shows an example of the inward radial ve- 
locity (-dr/dt) against particle radius at 5 AU in disc with 
a density of p(r) 



10 _11 (r/5AU)" 2 - 75 g cm -3 and a tem- 



perature of T(r) = 125(r/5AU)~"' J K. The radial velocity 
is clearly strongly dependent on the particle size, and can 
reach values of ~ 10 4 cm s _1 . 
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Figure 1. Radial drift velocity against planetesimal size at 5 AU 
in a disc with a density and temperature, at 5 AU, of 10 — 11 g 
cm -3 and 125 K respectively. 



The inward migration of dust particles, however, only 
occurs if the the pressure in the disc decreases mono- 
tonically with increasing radius. If there are any re- 
gions of the disc where the gas pressure increases with 
increasing radius, the gas velocity will become super- 
Keplerian, and the drag force will cause dust particles 
to move outwards towards the pressure maxima. This 
could occur in the presence of self-gravitating spiral s truc- 
tures (lHaghighipour fc Boss! 120031; iRice et all |2004) . in 
the presence of vortices iKlahTfcBod^nlielm^T2^03f) . or 
if an embedded planet has opened a gap in the disc 
JPaardekooper fc Mellemal 12004 12006ft . At the outer edge 
of a gap opened by a planet, the drag force will cause dust 
particles to migrate up the gap edge, away from the star. If, 
however, there is still gas accretion through the gap, which 
is expected for Jupiter mass companions around solar-like 
stars, then dust particles moving too slowly up the gap edge 
will be dragged by the gas into the gap. The gap edge will 
therefore act like a filter, allowing only particles of a certain 
size through into the inner disc. Even if the inner gas disc 
has not drained onto the central star, there could still be a 
substantial difference between the dust populations in the 
inner and outer disc and this could have a noticeable effect 
on the disc SED. 

To investigate this further, we need an approximation 
for the pressure gradient at the outer edge of a gap opened by 
an embedded planet. We therefore perform two-dimensional 
simulations of gaseous discs with embedded planets to de- 
termine the approximate disc structure. 



3 PLANET-DISC SIMULATIONS 

We simula te a planet embedded in a gaseous disc using the 
Zeus code iStone fc Normanll992T) . We use two-dimensional 
coordinates (r,cS) and a resolution of n r = 400, = 400. 
The computational domain extends, in code units, from 
rin = 0.5 to r out = 5 and we impose ouflow boundary con- 
ditions at both rin and r out . 

We assume that the disc is isothermal with a radially de- 
pendent sound speed given by c s (r) = 0.05r~ 1/2 . We adopt 
units in which G — 1 and the sum of the mass of the cen- 
tral star and embedded planet is M* + M p i = 1. Since the 



Figure 2. Disc surface density structure for the simulation in 
which the embedded planet has a mass of M p i = 10 — 3 . It is 
clear that the interaction between the planet and the disc gas has 
significantly depleted the surface density in the inner disc. 



disc thickness h/r is related to the sound speed through 
h/r ~ c s /Vk, we have a disc thickness of 0.05 at r = 1. We 
model angular momentum transport using a kinematic vis- 
cosity, v, that operates o nly on the azimuthal compon ent of 
the momentum equation iPapaloizou fc StanlevllT986ft . The 
initial disc surface density, E, is taken to have a radial de- 
pendence of oc r . Since M oc vYi(r) and since we 
would expect M to be constant we assume that, v(r) cx r. 
We normalise t he viscosity using the stan dard alpha formal- 
ism, v — ac 3 h dShakura fc Sunvaevlll973ri . and assume that 
at r = 1, a = 10" 3 . 

The planet is located at r — 1 and we consider planet 
masses of M p i = 5 x 10~ 4 , M p i = 10~ 3 , and M p i = 5 x 10~ 3 . 
The simulations are evolved for ~ 1000 planetary orbital 
periods. Figure [5] shows the surface density structure of the 
simulation with M p i = 10 -3 and illustrates that the inner 
regions of the disc have been significantly depleted. Figure 
[3] shows the surface density profiles for the three different 
planet masses. Also shown (dashed line) is an emperical fit 
to the gap edge. The surface density is normalised such that 
the total disc mass within r — 100 would be 10 -2 in code 
units. It appears that the chosen function is a reasonable 
fit to the gap edge, especially where the surface density is 
steepest, which is the region that would have the largest 
pressure gradient and hence would be the region where the 
drag force would have the largest influence on the embedded 
dust particles. From Figure[3]a general form for the function 
representing the gap edge would be 



-'edge 



E D r Ml- /exp 



-(r- 



>) 2 



d 2 

S a P 



(7) 



where r gap is the inner radius of the gap edge, and d gap 
is the gap scalelength. From the fits in Figure the gap 
scalelength is somewhat greater for M v \ — 5 x 10 -3 than 
for the lower planet masses, but we assume in general that 
dga P = r gap /5. The gap depth itself depends on the value 
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of the coefficient /. A large / value results in a deeper gap 
with a wider edge, and would correspond to a more massive 
planet. If we use Figure as a guide, / ~ 0.6 would corre- 
spond to a planet mass of Afpi = 5 x 1(T 4 , while / ~ 0.9 
would correspond to a planet mass of M p i = 5x 10~ 3 . These 
simulations have all been performed in scale free code units. 
To convert to real units, we simply need to introduce length 
and mass scales. For the problem we are considering here 
we will generally assume a mass scale of IMq giving a total 
disc mass of 1O _2 M0 and planet masses of about 0.5, 1, and 
5 Jupiter masses. The length scale we introduce depends on 
the assumed semimajor axis of the planet. In the rest of this 
paper we will assume a length scale of 4 au, resulting in a 
semimajor axis of 4 au and, from Figure [3] an outer gap 
edge at ~ 6 au. 



Planet moss = 5x10* 

1.9x10" 5 (au/r)(l-0.61«Exp(-(r-1.05)Vo.2 2 )) 3 



4 DUST MIGRATION AT GAP EDGES 

To determine how gas drag influences dust particles at the 
edge of a gap opened by an embedded planet, we need 
to know the gas pressure gradient. Equation Q gives an 
approximate form for the surface density structure at the 
edge of a gap. The volume density is then, p = E/ft where 
h — CsV/Vk- The sound speed, c s — kT / p,m,H where 
p = 2.4 is the mean molecular weight, k is Boltzmann's 
constant and t he temperature, T = 280(AU/r) 1/2 M,/M Q 
l)Havashi ITosl with M t = IMq. The pressure is then 
P — p(? s and the pressure gradient at the gap edge can be de- 
termined exactly since we have a fully analytic form for the 
pressure as a function of radius. We should, however, note 
that since the discs we consider here have optically thin inner 
regions, the outer edge of the gap will intercept more stellar 
flux a nd will be hotter than if the inner disc were optic ally 
thick JPullemond fc Domini3l2004t IXiceson et al.ll2005l) . A 
higher temperature will produce a larger pressure gradient. 
This temperature increase will, however, decay very quickly 
within the gap edge, and so may not be significant at the 
location where the pressure gradient is largest. Even if the 
temperature is slightly enhanced, a steeper temperature gra- 
dient will act to reduce the effect of the increased tempera- 
ture. Overall we would not expect the effect of the increased 
temperature at the edge of the gap to be particularly signifi- 
cant. Once the pressure gradient is calculated, the technique 
described in S|Ucan then be used to determine the radial ve- 
locity of solid particles at any point on the edge of the gap. 

The radial velocity of the solid particles will vary con- 
siderably along the edge of the gap. This is illustrated in 
Figure 0] which shows the radial velocity against particle 
size at the edge of a gap in which, r gap = 4 AU, d gap = 0.8 
AU, and / = 0.8. In this case the radial velocity is posi- 
tive, as opposed to Figure Q in which the radial velocity is 
negative. The surface density in Figure [I] is also normalised 
such that in the absence of a gap the total disc mass within 
100 AU would be 0.05M Q . The different curves in Figure 0] 
correspond to different locations (i.e., differents values of r) 
on the gap edge. All the curves in Figure 01 have the same 
form as in Figure At a given location, the radial veloc- 
ity intially increases as the particle size increases until it 
reaches a peak radial velocity, after which the radial veloc- 
ity decreases with increasing particle size. As r increases, 
the curves initially move upwards and the radial velocity in- 
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Figure 3. Disc surface density profiles for the three planet masses 
considered in the gasdynamic simulations. As expected, the gap 
depth increases as the planet mass increases. The dotted line also 
shows an approximate analytic fit to the profile of the gap edge 
which can be used to determine how the presence of a planet 
influences the interaction between the disc gas and the embedded 
solid particles. 



creases for all particle sizes (solid line, small dash line, and 
dash-dot line) reaching a maximum when r = r gap + 0.2 AU. 
As r increases beyond this point (dash-dot-dot-dot line, and 
long dash line), the curves move to the right (i.e., the particle 
size at which the radial velocity is maximal increases), and 
the maximum radial velocity decreases. The radial velocity 
for all particles smaller than the particle with the largest ra- 
dial velocity therefore decreases as r increases. Beyond the 
peak of the gap edge, the gradient changes sign, the radial 
velocity changes sign, and solid particles will drift inwards, 
rather than outwards. 
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Figure 4. Radial velocity against particle size at various loca- 
tions on the edge of a gap in which the gap edge is described by 
equation J7J with / = 0.8. As r increases, the radial velocities ini- 
tially increase reaching a maximum when r = r sa , p + 0.2au. As r 
increases further, the radial velocity decreases for particles smaller 
than the particle with the largest radial velocity (a < 1 — 10 cm). 
For particles small enough to contribute to the SED (a < 0.1 cm), 
there is therefore a location where they have a maximum outward 
radial velocity. 



What we wish to determine is the maximum outward 
radial velocity for particles that would contribute to the 
SED (eg., a <~ 0.1 cm). This can be determined by simply 
varying r until these particles achieve their maximum veloc- 
ity. In Figure 2] this would correspond to the dash-dot line 
(r — r gap + 0.2 AU). The two panels in Figure |S] show the 
maximum outward radial velocity as a function of particle 
size for two different surface densities and for three different 
/ values. In the upper panel the total disc mass within 100 
AU would be O.O5M0, while in the lower panel the total 
mass would be O.OIM©. In both figures r gap = 4 AU and 
dgap = 0.8 AU. The three curves in each figure correspond 
to / = 0.6, / = 0.8, and / = 0.9. Based on the simulation 
results shown in Figure |HJ these could correspond to planet 
masses of 0.5 Mj up , 1 Mj up , and 5 Mj up . 

The curves in Figure show that as /, or planet mass, 
increases, the peak radial velocity moves to smaller particle 
sizes. All particles smaller than the peak particle size then 
have increased radial velocities. The radial velocities in Fig- 
ure are, however, measured in the frame in which the gas 
radial velocity is zero. If there is no accretion (i.e. the radial 
gas velocity is zero relative to the central star) , all particles 
would be drifting up the gap edge, away from the central 
star. If there was some gas accretion, then particles with ra- 
dial velocities smaller than the radial velocity of the gas will 
be dragged into the gap. 

It is somewhat hard to quantify the radial gas velocity. 
The presence of a planet can significantly influence t he ve- 
locity structure of the gas in the vicinity of the planet dKlevI 
Il999l : iBate et, afll2003fl . We are, however, interested in the 
radial velocity of the gas as i t accretes t hrough the outer 
edge of the gap. Simulations (lKlevlll999t IBate et alJl2003t) 
show that in general the planet does not influence this ma- 
terial significantly, at least for planets with m asses below a 
few Jupiter masses IlLubow fc D 'Angelol2006Tl . The material 
that is most influenced by the planet is material that has al- 
ready accreted into the gap. Most of this material also ends 
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Figure 5. Radial velocity against particle size at the edge of a 
gap opened by an embedded planet. The two panels correspond 
to surface densities that differ by a factor of 5. The three curves 
in each panel correspond to different assumed gap depths, the 
deeper gaps assumed to result from the presence of a more massive 
planet. In both panels the peak radial velocity moves to smaller 
particles sizes as the gap depth increases resulting in an increase 
in the radial velocity of all particles smaller than the peak par- 
ticle size. The gap edge will therefore act to filter out the larger 
particles that have radial velocities larger than the inward radial 
velocity of the gas, allowing only the smaller particles into the 
gap. 

up within the planet's hill sphere, ultimately accreting onto 
the planet. The material that we are considering here is ma- 
terial that accretes through the gap into the inner disc, and 
it appears that this is largely unaffected by the planet, at 
least in the vicinity of the gap edge. The gas radial velocity 
through the edge of the gap will therefore be determined by 
the disc viscosity, and is related to the mass accretion rate 
through M = -2TrRT,v r . 

For an accretion rate of 10 -9 MQyr -1 the radial gas 
velocity in the absence of a gap would be ~ 1 cm s _1 for the 
higher of the surface densities (upper panel) in Figure|^|and 
5 cm s _1 for the lower of the surface densities (lower panel) 
in Figure |S] If the mass transfer through the gap edge is 
unaffected by the planet, the radial gas velocity should then 
increase as the surface density decreases. When / = 0.6, 
the surface density when the radial velocity of the parti- 
cles maximises is ~ 7 times lower than at the peak of the 
gap edge. The radial gas velocity should then be ~ 7 times 
higher, which for the upper panel of Figure |S] and for the 
mass transfer rate assumed above, would give v r ~ 7 cm s _1 . 
This means that only particles smaller than ~ 10 microns 
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will be swept into the gap by the radial motion of the gas 
since these particles have outward radial velocities less than 
~ 7 cm s _1 . For / = 0.9, the surface density changes by a 
factor of 860 and, if the same applies, then the radial gas ve- 
locity would then be ~ 860 cm s^ 1 . The maximum particle 
size that can populate the inner disc is then ~ 5 microns, 
comparable to that obtained when / = 0.6. An equivalent 
result is obtained if we consider the lower panel in Figure 
The surface density being 5 times smaller than in the upper 
panel results in the particles' radial velocities also increasing 
by a factor of 5. For a given mass accretion rate, however, 
the radial gas velocity also increases by a factor of 5 and 
hence the size of particles that would be swept into the gap 
is roughly indepedent of the assumed surface density. 

On the other hand, numer ical simulations by 
iLubow. Seibert fc Artvmowiczl <ll999T) suggest that the mass 
accretion rate through a gap actually decreases with in- 
creasing planet mass, at least for planets with masses in 
excess of 1 Jupiter mass. In this case the gravitational 
torques from the planet ac t to inhibit flow through the gap 
llLubow fc D'AngelolEoO^ . The rate of decrease with mass 
appears to be exponential, with the mass accretion rate de- 
creasing by almost an order of magnitude as the planet mass 
increases from 1 to 6 Jupiter masses. If this is indeed the 
case, then the size of particles that will be swept into the 
gap by gas accretion should decrease with increasing planet 
mass. Again considering the upper panel of Figure |S] if the 
three curves correspond to planet masses of ~ 0.5Mj up , 
~ lMj U p, and ~ 5Mj up then we may expect the gas radial 
velocity to decrease by an order of magnitude as / varies 
from 0.6 to 0.9. The maximum particle size swept into the 
gap would then decrease by an order of magnitude. For the 
gas radial velocities considered above, we might expect the 
maximum particle size to decrease from ~ 10 microns to a 
few tenths of a micron as / increases from 0.6 to 0.9 corre- 
sponding to the planet mass increasing from 0.5 Mj up to 5 
M Jup . 

Since the above process results in only the smallest dust 
grains being swept into the gap, we would expect the dust 
to gas mass ratio to be different in the inner disc, compared 
to the outer disc. It is, however, difficult to make an inde- 
pendent measurement of the gas and dust masses, and hence 
determining the different dust to gas mass ratios is probably 
not possible. The presence of a silicate feature in protostel- 
lar disc spectral energy distributions (SEDs) does, however, 
depend on the grain size distribution. The presence of large 
grains tends to wash out the silicate feature. As we discuss 
in the next section, SED modelling of those systems with 
near infrared deficits suggest that the grain size distribution 
in the inner disc differs from that in the outer disc in a way 
that is consistent with dust filtration at a gap edge. 



5 DISCUSSION 

Although the suggestion that the observed near-IR deficits 
in some TTauri SEDs could result from an embedded 
Jupiter-like planet clearing the inner disc is exciting, there 
is, as yet, no clear evidence that this is indeed the case. 
Other non-planetary possibilities, for example rapid grain 
growth, could also result in optically thin regions that may 
produce SEDs largely consistent with that observed. 



A sufficiently massive planet does, however, have a real 
effect on the structure of the gaseous disc which may influ- 
ence the distribution of solid particles. This effect could, 
in turn, influence the nature of the resulting SED. A ~ 
Jupiter mass planet should open a gap in the disc, the in- 
ner region of which could clear to form an inner hole. The 
outer edge of this gap (or hole) will have a pressure that 
increases with increasing radius, increasing the net outward 
force on the gas and resulting in super-Keplerian gas veloc- 
ities. The resulting drag force on the solid particles then re- 
sults in these particles drifting up the gap edge, away from 
the central star. If, however, there is still some gas accre- 
tion int o the gap - which is expected for ~ Jupiter-mass 
planets llLubow. Seibert fc Artvmowiczlll999l) - solid parti- 
cles with radial velocities smaller than the radial gas velocity 
will populate the inner region of the disc. For disc properties 
appropriate for TTauri stars, the radial velocity increases 
with particle size for particles smaller than ~ 1 cm. This 
means that only the smallest particles will have outward ra- 
dial velocities small enough to be dragged into the gap by 
the accreting gas. 

Gas accretion through the gap edge therefore acts to 
filter solid particles, with only the smallest particles making 
it into the gap. The exact size range that is influenced by 
the gas depends on the properties of the gas disc and on the 
radial gas velocity. Assuming a functional form for the gap 
edge (motivated by two-dimensional, hydrodynamic simula- 
tions of embedded Jupiter-like planets) we show that if the 
inward gas velocity is ~ 1 — 10 cm s -1 , only ~ micron and 
smaller grains will be able to populate the inner disc. If the 
gas accretion rate is modified by sufficiently massive planets, 
as su ggested by simulations llLubow. Seibert fc Artvmowiczl 
1999), the size of particles that may populate the inner disc 
should decrease with increasing planet mass. Systems, such 
as CoKu Tau/4, that appear to have almost no material in 
the inner disc may indicate the presence of a massive com- 
panion while those, such as GM Aur and TW Hya, that 
appear to still have some material in the inner disc, may 
have a companion with a mass that allows some material to 
continually replenish the inner disc. 

What is attractive about this possibility is that SED 
modelling suggests that system with near-IR deficits do in- 
deed have two populations of dust particles. This spatial 
segregati on of dust within di scs has been inferred in the TW 
Hya disc dCalvet et al.ll20o3) , and radiation transfer models 
of new Spitzer Space Telescope IRS spectra of the CoKu 
Tau 4, GM Aur and DM Tau discs also require inner holes 
with small grains in the inner disc, and a larger grain popu- 
lation in the outer disc jD'Alessio et al.ll2005l : ICalvet et alJ 
2005). The silicate feature present in the IRS spectra provide 
the clues that the dust grains in the inner disc are smaller 
than in the outer disc. As more data comes in from Spitzer, 
the number of disc systems showing evidence for inner holes 
is increasing (lLada et all 120061: iMuzerolle et afll2006h . To 
date, most candidates have IRAC and MIPS photometry, 
so only the presence of inner disc clearing can be inferred. 
Follow up observations with IRS will enable us to determine 
whether the spatial segregation of dust sizes is a common 
feature and thereby further test our dynamical models for 
disc-planet interactions and dust filtration. 

The SED modelling of systems like GM Aur and DM 
Tau not only require a population of small grains in the inner 
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disc, but also require that the mass, or surface density, is also 
low enough for the inner disc to be optically thin. Accretion 
onto a planet can reduc e the amount of mass rea ching the 
inner disc by up to 90% llLubow fc D'Angeloll200rJ) . but this 
is probably insufficient to produce an optically thin inner 
disc. Typically the optical depth will need to be reduced by 
a factor of ~ 10 4 , a factor of 1000 more than can be produced 
by the planet alone. If solid particles are filtered at the gap 
edge, the amount of mass reaching the inner disc will depend 
on the size distribution of the solid particles. Interstellar 
medium (ISM) grains are generally assumed to have sizes 
betwe en 0.005u-m and lpm and a size dis tribution of n(r) oc 
r~ 3 ' 5 llMathis. Rumpl fc Nordsieckll977l) . Grain growth can 
however change this distribution quite significantly. If there 
is constant rep lenishment of the ISM grains, this ca n steepen 
to n(r) oc r -4 iMizuno. Markiewicz fc Voelkll988l) . in which 
case most of the mass will be in small grains. If, however, 
there is no source of ISM grains, the final size distribution 
can be much flatter, approaching n(r) oc r~ 2 . A relatively 



shallow size distribution (e.g.. 



will already have a 



reduced opacity, independent of the presence of a planet, 
but, as shown below (see also Figure[7|l, a further reduction 
is needed to match the observed SEDs. Such a reduction can 
be obtained through the filtration process described here. 

The amount of material that can reach the inner disc 
will depend on the maximum size that can accrete through 
the gap edge, and on the size distribution of the particles 
in the outer disc. Figure [(J shows the ratio of the mass of 
particles smaller than lu.m (solid line) and 10 pm (dashed 
line) to the total mass of solid particles, plotted against 
the exponent of the power-law size distribution. In this fig- 
ure we assume a minimum size of 0.005 pm and a maxi- 
mum size of 1 mm. For steep size distributions most of the 
mass is in small grains, and even if the maximum size of 
filtered particles is lpm, the inner disc would remain op- 
tically thick. It is also likely that in this case, grain coag- 
ulation would soon produce large grains in the inner disc 
JPullemond fc Dominikll2005l) . From SED modelling, how- 
ever, we know that systems like GM Aur, TW Hya and DM 
Tau must have a size distribution that extend to larger sizes 
in the outer disc and, since these systems are reasonably 
evolved, there is unlikely to be a significant new source of 
ISM grains. We might, therefore, expect the size distribution 
in the o uter disc to be flatter than the ISM grain size distri- 
bution iMizuno. Markiewicz fc Voelklll988f) . It has already 
been shown that in some systems, flatter grain size distri- 
bution (e.g., n(r) oc r~ 2 ' 5 ) provide better fits to the SEDs 
than ISM size distributions iD'Alessio. Calvet fc Hartmannl 
l200ll) . Figure [5] shows that a size distribution of n(r) oc r -3 
and a maximum filtered size of lpm, would reduce the mass 
of solid material accreting through the gap edge by a fac- 
tor of 10 3 . If the planet accretes 90% of this material, the 
amount of solid material reaching the inner disc would be 
reduced by a factor of 10 4 , sufficient to produce an optically 
thin inner disc. 

This process would also produce a very different dust- 
to-gas ratio in the inner disc, compared with the outer disc. 
That systems like GM Aur and DM Tau have optically thin 
inner discs but are still accreting at M 10 -9 MQyr -1 sug- 
gests that there is still a substantial amount of gas in the 
inner disc. The optically thin inner disc requires a reduction 
in solid particles of order 10 4 , while the ongoing accretion 
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Figure 6. Ratio of the mass of particles with sizes below lpm 
(solid line) and 10 pm (dashed line) to the total mass of solid 
particles, plotted against the exponent of the power-law size dis- 
tribution. If the maximum size of particles filtered at the gap edge 
is 1pm, and if the power-law exponent is —3, the mass of material 
accreting into the gap would be reduced by a factor of 10 3 . If the 
planet accretes 90% of this mass, the amount of material reaching 
the inner disc would be reduced by 10 4 , sufficient to produce an 
optically thin inner disc. 



requires a far smaller reduction in the amount of gas in the 
inner disc. This is at least qualitatively consistent with the 
suggestion that the solid particles are filtered at the outer 
edge of the gap which, together with the subsequent accre- 
tion onto the planet reduces the mass of solid particles by 
a factor of ~ 10 4 , while at least 10% of the gas is a llowed 
through into the inner disc llLubow fc D 'AngelcfeOOd) . Such 
a reduction in the solid particle surface density in the inner 
disc also means that sub sequent grain coagulation is p roba- 
bly not very important. iDullemond fc Dominikl (2005). us- 
ing a relatively simple one-particle model iSafronovl ll969). 
show that at ~ 1 au grain coagulation can produce ~ 500pm 
particles within a few hundred years for reasonable surface 
densities. If the mass that can be swept up by a growing par- 
ticle is then reduced by a factor of 10 4 , one might assume 
that the maximum size that this particle could achieve is 
then reduced by a factor of 20 (i.e., a maximum size of ~ 25 
microns). This, however, is a significant overestimate since 
the amount of mass swept up by a growing particle depends 
strongly on the size of the particle, or its collision cross- 
section. In the simple one-particle model, there is almost no 
grain growth if the surface density is reduced by a factor 
of 10 4 from an initial density that would be appropriate for 
the inner regions of a TTauri disc. Although this is a fairly 
simple model, it at least suggests that if the amount of mass 
reaching the inner disc is reduced significantly, coagulation 
should not play an important role in the subsequent evolu- 
tion of these grains. 

Figurc|7] illustrates how the filtration process discussed 
above may influence the disc SEDs. The SED models in 
Figure |7| are computed assuming a 300 au flared disc of to- 
tal mass O.OIMq that is passively heated by a star with 
Tstar = 4000 K, and i? s tar = Rq- The disc surface density 
and scaleheight have radial profiles of £ oc r _1 and h oc r 1 ' 25 , 
and the scale height is normalised such that h = 18 au at 
r = 100 au. The temperature is determined self-consistently 
by the Monte-Carlo radiation transfer calculation. The solid 
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Figure 7. SED models that illustrate how filtration of dust par- 
ticles at the outer edge of a gap can influence the resulting SED. 
These models are computed assuming a 300 au flared disc of total 
mass O.OIMq that is passively heated by a star with T s tar = 4000 
K, and i? s tars = Rq- The solid curve is the SED resulting from 
a dust model with a grain size distribution of n(r) oc r -3 with 
minimum and maximum grain sizes of 0.005 um and f mm respec- 
tively. The other two models have the same grain population in 
the outer disc, beyond 5 au, but have fSM-like grains in the inner 
disc with depletion factors of f0 4 (dashed line) and 10 6 (dotted 
line). Notice the near-fR deficit and prominent silicate features 
at f0|am and 20um for the SED models with ISM-like grains and 
inner disc depictions, consistent with observations of a number of 
T Tauri systems. 

curve is the SED resulting from a dust model with a grain 
size distribution of n(r) oc r -3 with minimum and max- 
imum grain sizes of 0.005(xm and 1 mm respectively (see 
IWood et alJ <2002ft . Table 1 and Figure 3). The other two 
models simulate the dust filtration discussed above. These 
models have the large grain model in the outer disc (be- 
yond 5 au) and ISM-like grains (maximum size lum) inside 
5 au. The dashed curve has a depletion factor of 10 4 for the 
dust density, while the dotted curve has a depletion factor of 
10 6 , still within the possibilties discussed above. Notice the 
near-IR deficit and prominent silicate features at 10u.m and 
20[im for the SED models with ISM-like grains and inner 
disc depletions. The redistribution of IR excess from short 
to long wavelengths for disc s with optically th in interiors is 
also apparent as discussed in lRice et alJ <l2003t) . These SEDs 
are not models for any particular system, but demonstrate 
the observable signatures of the filtration process we have 
described. In these models, the dust density in the inner disc 
is significantly reduced by filtration at the gap edge and by 
accretion onto the planet, while the gas is only influenced 
by accretion onto the planet. The gas density in the inner 
disc is therefore reduced by a much smaller amount than the 
dust density. This extremely small dust-to-gas ratio in the 
inner disc produces an "opacity gap" , while still allowing gas 
accretion onto the star, in line with observations of many T 
Tauri stars that exhibit near-IR deficits in their SEDs. 

It has, however, also been suggested 
feisner. Chiang fc Hillenbrand! lioofj that "opacity gaps" 
and the presence of small grains in the inner disc could 
actually be evidence for grain growth and subsequent rapid 



inward migration, rather than inferring the presence of a 
gap opening planet. This conclusion was partly because 
Jupiter-like planets in relatively massive discs should 
migrate into the central star wit hin the disc lifetime 
feisner. Chiang fc Hillenbrand! I20061) . Although this may 
indeed be the case, the viscous timescale at 5 au can be 
~ 10 6 years for reasonable viscosity values. The observation 
of such a system is therefore still possible, even if the planet 
ultimately doesn't survive the migration process. If near-IR 
deficits are due to grain growth, there must also be some 
means of sharply truncating the grain size distribution at 
~ lum in the inner disc. It is possible that we are seeing 
a population of planets that form early, and potentially do 
not actually survive the migration process. 



6 CONCLUSION 

Although we cannot definitively conclude that near-IR 
deficits in the SEDs of some TTauri stars are a consequence 
of an embedded Jupiter-like planet, that the filtering effect 
of such a planet - producing a small grain population in 
the inner disc that differs from the population in the outer 
disc - is consistent with SED modelling of these systems is 
suggestive. The SED modelling also requires that the inner 
disc be optically thin, while ongoing gas accretion requires 
that the inner disc still have a substantial amount of gas. 
If the grain size distribution i s somewhat flatter than the 
stand ard ISM size distribution dMathis. Rumpl fe Nordsieckl 
1977), filtration at the outer gap edge can significantly re- 
duce the amount of solid material reaching the inner disc, 
while still allowing sufficient gas through the gap. Although 
we don't actually know the grain size distribution in the 
outer discs of those systems that show inner holes, they 
do show evidence for grain growth. In some cases flatter 
size distributions do provide better fits to TTauri SEDs 
iD'Alessio. Calvet fc Hartmannl200Jl . although this has yet 
to be tested on the systems we consider here. That filtration 
of solid particles at the outer egde of a gap opened by an 
embedded Jupiter-like planet can produce not only a popu- 
lation of small grains in the inner disc, but can potentially 
also make this region optically thin, makes this an attractive 
possibility. 

The filtration process should also lead to an enhance- 
ment of larger particles at the peak of the gap edge. Parti- 
cles exterior to the gap will continue to drift inwards until 
they reach the gap edge, at which point the change in the 
gas pressure gradient will prevent the larger particles from 
drifting any further. If the surface density of such particles 
becomes sufficiently high, further planet formation could oc- 
cur at this location. Since 'Hot' Jupiters are thought to form 
at modest radii and then migrate inwards via Type II mi- 
gration to their present locations, one may expect in some 
cases to find additional lower mass planets in 2:1 resonances 
with these 'Hot' Jupiters. The presence of even terrestrial 
mass planets in resonant orbits with 'Hot' Jupite rs could be 
detected using transit timing jAgol et al.ll2005l) . Although 
such a detection would be consistent with the analysis pre- 
sented here, it has been suggested that this could also occur 
if low-mass, fast migrating planets (> 0.1 Mg,) were r eso- 
nantly trapped by gap opening planets (iThomme IHH). In 
the latter case, however, one may expect the trapped plan- 
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ets to subsequently grow into gas giants jThommesl [2005 ) . 
while in the former case the planetesimal trapping could 
occur later when there may only be sufficient time for the 
newly forming planet to reach a terrestrial-like mass. 
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